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ABSTRACT 

The breakout of a supernova shock wave through the progenitor star's outer envelope is ex- 
pected to appear as an X-ray flash. However, if the supernova explodes inside an optically- 
thick wind, the breakout flash is delayed. We present a simple model for estimating the con- 
ditions at shock breakout in a wind based on the general observable quantities in the X-ray 
flash lightcurve: the total energy Ex, and the diffusion time after the peak, ^(jiff ■ We base the 
derivation on the self-similar solution for the forward-reverse shock structure expected for 
an ejecta plowing through a pre-existing wind at large distances from the progenitor's sur- 
face. We find simple quantitative relations for the shock radius and velocity at breakout. By 
relating the ejecta density profile to the pre-explosion structure of the progenitor, the model 
can also be extended to constrain the combination of explosion energy and ejecta mass. For 
the observed case of XRO08109/SN2008D, our model provides reasonable constraints on the 
breakout radius, explosion energy, and ejecta mass, and predicts a high shock velocity which 
naturally accounts for the observed non-thermal spectrum. 

Key words: shock waves-supemovae: general-supernovae: individual: SN 2008D stars: 
winds, outflow 



1 INTRODUCTION 

As a supernova shock wave propagates through the progenitor star, 
it eventually emerges through the outer envelope region which has 
a low optical depth. At this point the shock "breaks out" of the 
star, initiating the first electromagnetic signal of the explosion. 
The breakout results in a flash of radiation, as the internal energy 
deposited by the shock diffuses out of the shocked region on a 
timescale comparable to the dynamical time of the shock. Follow- 
ing breakout the ejected envelope expands so that the photoshpere 
recedes into the ejecta. Any remaining internal energy, which has 
not been converted during adiabatic expansion, is then gradually 
radiated to power the supernova lightcurve. 

It has long been suspected that the shock breakout from a star 
would appear as an X-ray flash ( |Colgate|1974[ jKlein & Chevalier| 
|1978[ >, followed by a UV/optical transient corresponding to emis- 
sion from the outer layers of the ejecta ([FaK|[T978]l. The interest 
in prompt signals from supemovae increased recently due to the 
new capabilities of modem searches for transients, such as Pan- 
STARR^and PTlj^as well as the planned LSSlj^and WFIRSlQ 
In fact, several supernovae lightcurves were observed early after the 
explosion both in the UV and the optical bands l |Gezari et al.|2008] 
[20T0||Soderberg et al.|2008^ . In two cases, GRB060218/SN2006aj 

^ http://pan-stan's. ifa.hawaii.edu 

^ http://www.astro.caltech.edu/ptf/ 

^ http://www.lsst.org/lsst 

^ http://wfirst.gsfc.nasa.gov/ 



(|Campana et al.|2006| l and XRO08109/SN2008D jSoderberg et"d^ 
'2008V X-ray detectors on board the Swift satellite have captured 
a luminous X-ray outburst, which were later followed by observa- 
tions of supemovae in the longer wavelengths. Both events may 
very well have been the signature of shock breakout. 

These new observations motivated detailed theoretical model- 
ing of the early emission in supemovae. The models focused on the 
early UV/optical lightcurve that follows breakout, on time scales 
of or der a day jChevalier & Fransson|2008[ [Rabinak & Waxman| 
|2010[[Nakar & S ari'2010V for which the dynamics are somewhat 
simpler (see also Che valier] ( |1992[ l). Modeling the X-ray flash as- 
sociated with breakout is more complicated, since it must include 
the properties and structure of the radiation-mediated shock (RMS) 
as it propagates through the sharply decreasing density of the outer 
envelope, where the shock width becomes non-negligible relative 
to the distance to the edge of the star. jKatz et al.| ( (2010^ have shown 
that if the shock velocity normalized by c, Ps 55 0.07, free-free 
emission will be unable to produce a large enough number den- 
sity of photons to establish a blackbody spectrum. As a result, the 
average photon energy in the shocked region will be significantly 
larger than the equilibrium blackbody photon temperature with the 
same energy density. Hence, any analysis of the breakout lightcurve 
based on the assumption of local thermodynamic equilibrium and a 
black body spectmm will not reproduce the photon flux and spec- 
trum emitted from the shock region during breakout. 

The shock breakout lightcurve could be significantly influ- 
enced by any circumstellar material (CSM) surrounding the explod- 
ing star. In fact, the aforementioned two candidates for an X-ray 
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flash of a shock breakout origin are associated with a type Ibc su- 
pernova, in which the progenitor is beheved to be a Wolf-Rayet 
star that has undergone significant mass loss prior to the explosion. 
General, order of magnitude estimates for the energy enclosed in 
shock breakout do indeed favor a wind environment over a bare 
star in both cases ( fWaxman et al.|2007[pievaUer & Fransson|2008| 
|Katz et al.|2010^ . If the supernova takes place in an optically thick 
wind, the breakout is delayed from the observer's point of view, 
and the X-ray flash emerges only when the shock width in the wind 
becomes comparable to its distance from the CSM photosphere. 
Even a dynamically unimportant wind (which does not influence 
the emission of the early UV/optical lightcurve), will completely 
change the properties of the X-ray flash relative to a bare star by 
virtue of its optical depth. 

Motivated by these considerations, we present here a simple 
analysis that relates the observable properties of the X-ray flash 
to the physical parameters of shock breakout in an optically thick 
wind. While a self-consistent treatment of the shock structure, light 
travel time and frequency dependence is required to reproduce the 
lightcurve of the flash, the integral quantities should be weakly de- 
pendent of these details. We apply the known self similar solution 
of a supernova ejecta moving in to a pre-existing wind material with 
a density profile of p oc r^^ ( Chevalier||l982 and demonstrate 
how it can be used to relate the key observable features of the X- 
ray flash, the total energy Ex and the diffusion time scale t^is to 
the parameters of the shock breakout. Most notably, we can obtain 
a quantitative estimate of the breakout radius and the associated 
shock velocity. Our model expands upon the order of magnitude 
estimate presented by |Ofek et al.| ( |2010) l regarding PTF09uj, which 
possibly exhibited a long (~ 1 week) UV flash due to breakout in a 
very dense wind. While our model assumes a non relativistic shock, 
it can be generalized to relativistic shocks and more complex wind 
structures, and also compared to the specific model suggested by 
[LilpOO?) for SN2006aj (note, however, that this source may have 
been significantly asymmetric, (Waxman et al.|2007 1). 

The outline of the paper is as follows. In ^ we review the self 
similar solution for an ejecta expanding into a wind, and discuss 
the applicability of this physical setting to the problem of shock 
breakout. In S|3]we present how the self similar solution can be used 
to relate the observed properties of the X-ray flash to the conditions 
at breakout. In ij4]we demonstrate that based on the initial profile 
of the ejecta ( |Matzner & McKee|1999| l, the observable quantities 
can be used to place constraints on the parameters of the underlying 
explosion. In [|5]we consider the particular example of SN2008D, 
which is a borderline case in terms of the applicability of our model. 
Finally, !|6]summarizes our main conclusions. 



2 THE FORWARD-REVERSE SHOCK STRUCTURE 

We consider a spherically symmetric supernova explosion in a 
nearly stationary CSM, produced by a wind from the progenitor 
star. The wind is assumed to be much slower than the ejecta. For a 
steady wind with a mass loss rate M and a velocity Vw , the density 
profile outside of the star is. 
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Following the explosion, the outgoing ejecta plows through 
the wind material and slows down, driving a forward shock through 
the wind and a reverse shock that propagates back into the ejecta. 
Once the shock had propagated far enough from the star's surface 



10 




Figure 1. The velocity (black), pressure (red) and density (blue) profiles of 
the self similar solution for supernova ejecta propagating through a wind 
with m = 10 and 7 = 4/3 (see text). The radial coordinate is normal- 
ized to the radius of the forward shock, R p , and the physical quantities are 
normalized to their post shock values at the forward shock front. 



so that the crossing time and the size of the star can be neglected, 
one may assume an asymptotic time-dependent density profile for 
the outer part of the supernova ejecta of the form. 



(2) 



where g and m are constants which depend on the initial conditions 
of the progenitor and the explosion. If both the ejecta and the wind 
material are initially cold, a self-similar solution can be found to 
describe for the structures of both the forward shock l iParker|1963[ > 
and the reverse shock ( |Chevalier|1982^ . At any time, the flow con- 
sists of a forward shock region of the swept up wind material, and 
a reverse shock in a leading region of the ejecta. At the bound- 
ary between the two regions, velocity and pressure are continuous 
but density is not. Denoting the position of the forward shock by 
Rf, the reverse shock by Rr and the contact radius by Rc, the ra- 
tios Rr/Rf and Rc / Rf are time-independent. The solutions are 
found by integrating dimensionless functions of similarity variable, 
77 — r^^'^t^^, where A = (n — 2)/ (n — 3). The forward and reverse 
shock fronts serve as boundary conditions. Integration is carried out 
up to Rc from both ends, and the solutions are joined by matching 
the velocities and pressures on both sides of the contact surface. 

The self-similar solution describes the normalized profiles of 
the density, velocity and pressure profiles as a function of 77, and 
depends on the values of m and the adiabatic index of the shocked 
materials, 7. An example of such a profile for m = 10and7 = 4/3 
is presented in Figure [T| where all quantities are normalized to 
their values just behind the forward shock. The gradual divergence 
of density at the discontinuity (in both of the forward and reverse 
shock regions) is generic to the assumed wind profile ('Cheva-| 
|lier| |1982|, for which the solutions lead to a vanishing speed of 
sound on both sides of the contact radius. In reality, the divergence 
will be avoided since the density discontinuity is Rayleigh-Taylor 
unstable and hence the contact region will obtain a finite width. 

Our choice of 7 = 4/3 is, of course, motivated by our assump- 
tion that the shocks involved are radiation dominated; the internal 
energy density in the shocked region is therefore e = 3P, where P 



is the pressure (see I Chevalier 1983 1 for comparison, with 7 = 4/3 
motivated by cosmic-ray dominated shocks). By integrating over 
the pressure in the solution we can find the total energy in radia- 
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tion. The pressure behind the forward shock, Pf, is determined by 
the strong shock condition, 



Pf = 



7 



(3) 



where Rf is the velocity of the forward shock. Denoting the rela- 
tive width of the shocked region by a:: = (Rf — Rr) /Rf we may 
express the total internal energy in the shocked region as: 
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Xe^i- j^RfxRp. 



The dimensionless quantity xe is defined by 



XE = 



(4) 



(5) 



4-kR%xPf 

and is determined by the details of the post-shock profile. For the 
particular case of m = 10 and 7 = 4/3 depicted in Figure [T] 
X = 0.146, XE = 0.9. 

The total width x is dominated by the forward shock region, 
which has a width of xf = {Rf — Rc)/Rf = 0.135, while the 
width of the reverse shock region is only xr = {Rc — Rr) /Rf ~ 
0.011. However, the density in the reverse shock is much larger, and 
we find that Mr/Mf = 2.23, where Mr and AJf are the masses 
enclosed in the reverse and forward shock regions, respectively. 
By construction, Mf = 4:TtKRf- The resulting {AIr + Mf) ~ 
3.2A/f implies that order of magnitude estimates identifying the 
total energy of the flash with MfR^f/"^ (see, e.g., |Katz et al.| 
( |2010j^ ) should be corrected. While it is clear that the mass of the 
ejecta affected by the wind must be comparable to that of the swept- 
up wind, ignoring the mass enclosed in the reverse shock leads to a 
significant overestimate of the shock velocity. 

The density profile in the self similar solution can be inte- 
grated to find the normalized column density in the forward and 
reverse shock regions, which can then be used to estimate the opti- 
cal depth and the diffusion time scales. First, we define. 
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rKRp^'iRF- Rc) ' 7KRf{Rc-RR)' 
The results for the self similar solution depicted in Figure [T] are 
tf — 1.3 and tr ~ 4.2. The total optical depth of the shock 
region is therefore 



Ts = K, {tfxfRf + trXrRf) 7KRf^ — 1.55k,K Rf 



(7) 



where we have assumed a uniform, frequency independent opacity, 
K. The total diffusion time scale in the shocked region, tus is then. 



tos = Ts- 



xRf 



0.23 



(8) 
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Finally, the physical scales of the full solution are determined 
by the dimensional parameters, g and K. In particular, the velocity 
of the forward shock is related to its position by, 

l/(m-3) 



Rf = 
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(9) 
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where j4f is another constant that is extractable from the self- 
similar solution I Chevalier|1982 1; for m = 10 and 7 = 4/3 we 
get ^ 0.07. 

We note that the numerical values of x, xe, t'f and tr are 
very weakly dependent on m. In the range of m = 8-12, we find 
that these values vary by no more than ±20% with respect to their 
values for m = 10, which we use for reference below. The value of 
Af is not as well constrained, and changes within a factor of 2 in 
this range. However, as shown in §|4] the final result is not sensitive 
to this uncertainty. 



3 A SIMPLE PARAMETRIZATION OF THE BREAKOUT 
CONDITIONS 

We now use the results of the previous section to estimate the shock 
breakout features. We assume that prior to shock breakout the flow 
is approximately adiabatic in the shocked regions, and that leakage 
of radiation into the cold parts of the wind and ejecta is negligible. 
We discuss the validity of these assumptions below. 

We take the total energy observed in the breakout flash. Ex , to 
be comparable with the energy carried in radiation in the shocked 
region, E^ad- Similarly the time scale characterizing the post max- 
imum part of the lightcurve is determined by diffusion from the 
shock front through the unshocked wind, tdift - These associations 
are valid because the diffusion time of radiation in the shocked re- 
gion is shorter than the time scale for diffusion from the shock to the 
unshocked wind's photosphere. The optical depth from the shock 
front in the wind material is r = kK/ Rp, and so the time scale 
for diffusion from the shock front is 



Rf kK 



(10) 



or about four times longer than the timescale for diffusion in the 
shocked region, given by equation The opposite regime of 
tdiff tos would have been characterized by a longer, fainter 
lightcurve. 

Observationally determined values of Ex and t^ia and the 
self similar solutions provide the necessary relations for estimating 
the key properties of ejecta-wind configuration at breakout. First, 
using Equation \10\ the observed value of tdiff is sufficient to de- 
termine the value of the wind density constant, K: 



K 



= 1.5 X 10 — g cm , 

K, K0.2 



(11) 



where tdiff.a ~ (tdiff/10^ s) and K0.2 = (k/0.2 cm^ g~^) (the 
fiducial opacity normalization is appropriate for wind material that 
is dominated by helium and heavier elements.) The value of K 
for tdiff,3 = 1 corresponds to a wind mass loss rate of about 
3 X 10^'' Mq yr~^ for a wind velocity of 100 km s'^. 

The criterion for shock breakout is that the diffusion time be 
equal to dynamical time, tdift ~ tdyn = Rf/Rf- We use this 
relation to substitute the shock velocity from Equation |4]l for tdiff, 
yielding 
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Using the predetermined values of xe and x based on the self- 
similar solutions in [j2] the breakout radius can be expressed in 
terms of Ex and tdUf '■ 
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5.4 X 10' 



ExAltdiB;, 
K.0.2 



1/3 



0.146 



-1/3 



V0.9/ 



-1/3 



(13) 



where Ex, 47 = (-Ex/ 10*'^ ergs). We emphasize again that the 
specific values of x and x e found in the self-similar solutions are 
weakly dependent on the exact value of the ejecta density profile, 
771. Combined with the fact that the breakout radius depends only 
on the cubic roots of these numerical factors, we conclude that the 
the result of Equation \13\ is robust to uncertainties in the details 
of the ejecta profile. 

Using the inferred values of the breakout radius and the wind 
density coefficient, K, we can also express the masses of the wind 
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material and the ejecta included in the forward and reverse shock, 
respectively. The mass in the forward shock is: 

Mf _ 

Mo " 



5 1 X in-^/?^/^ f^/^ ,,-4/3/ ^ y 



(14) 



JZ^\-V3 /Xe\ -1/3 

.0.146/ V0.9/ 
and Mr is about 2.23 times larger. It is a useful consistency check 
that the mass of the ejecta involved in the shock formation is signif- 
icantly smaller than the total ejecta mass. An ejecta density profile 
of the form of equation (j2]l is only applicable to the edge of the 
progenitor envelope, where the initial density drops sharply with 
radius. 

An important consequence of the estimate for the breakout ra- 
dius is the corresponding shock velocity. The value of this velocity 
(normalized by c) is, /3s ~ O.IS-E"^^^ 47^diff'3- This relates the ob- 
served global features of the lightcurve and the spectrum, since for 
/3 > 0.07 the photon generation rate in the shocked region is in- 
sufficient to establish a Planck distribution within the shock width 
( |Katz et al.|2010^ . The emergent spectrum in that regime is dom- 
inated by harder IkeV) X-ray photons than a blackbody, since 
the same amount of energy is shared by fewer photons. If, on the 
other hand, Ps < 0.07, the flash should exhibit a blackbody spec- 
trum with a temperature of 
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2.64 X W'Elil^.t^,^ 
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/12 1/6 / 

,3 '^0.2 



1/6 /XB\l/6 



.0.1467 vo.gy ^' ^''^ 

assuming that the characteristic luminosity is L — 

A-kRIo^T^b ~ Ex/Um- 

Some cautionary remarks should be mentioned about our re- 
sults. First, a useful aspect of equation l |13[ l is that for the refer- 
ence values of Erad ~ 10'*^ergs and tdiff ~ lO^s, the light cross- 
ing time at the breakout radius is about 18% of the diffusion time. 
Therefore, the geometric time delay due to the light crossing time 
will have a secondary, but not negligible, effect on the earlier part 
of the flash lightcurve. Another caveat is related to our assumption 
that the kinetic energy of the shocked region (and that of the un- 
shocked ejecta) will not be converted into internal energy over the 
time scale of the breakout flash. We expect this approximation to 
be valid, since the shocked region would slow down significantly 
only after it sweeps up an additional wind mass comparable to its 
own, {Mr + Mf) = 3.2 x AtvKRbo (see Sj2](. Hence, it takes 
at least several dynamical times, or diffusion times, to convert ki- 
netic energy to internal energy. Moreover, the kinetic energy of the 
shocked material, Ekin ~ ^{Mr + Mf){Rf/7)^, is only about 
10% of its internal energy, and so the addition to the available en- 
ergy would generate at most a 3% correction in the value of _Rso in 
equation l |13^ . The main reservoir of kinetic energy lies in the yet 
unshocked ejecta, but this region takes even longer to slow down 
and is initially hidden behind larger optical depths. 

Another secondary effect we have neglected is the possible 
acceleration of gas ahead of the shock by radiation which escapes 
the shocked region. While an exact calculation of the extent of this 
effect requires a radiation-hydrodynamics simulation, it is straight- 
forward to demonstrate that relative to velocities > lO^cm 
significant acceleration of wind material at any given radius be- 
low Rbo cannot take place before the shock overtakes that radius. 
It is also noteworthy that the optical depth of the unshocked ejecta 
is quite large, so radiative effects should not increase the width of 
the reverse shock significantly prior to breakout. 

A more serious concern is to what extent the ejecta is in a 



cold, power law density density state prior to breakout. For com- 
pact progenitors with a radius i?* ~ lO^^cm, the reference values 
give Rbo/R* ~ 54. While this ratio is sufficiently large to justify 
neglecting the progenitor radius as a relevant scale in the solution, 
it is not obvious that the density profile should already settle at 
Rbo into the asymptotic form we assumed in equation |2|. Origi- 
nally, the supernova shock wave deposits in the outer envelope ~ 6 
times as much internal energy as kinetic energy. For spherical adi- 
abatic expansion, the internal energy declines as ~ r^^, so for 
{EradA7'tdiB,3) < 1 the ejecta can still be "lukewarm" rather than 
cold, implying that the ejecta density profile may still be develop- 
ing, and perhaps that the reverse shock may not be very strong. This 
might lead to some deviation of the gas dynamics from our simple 
model. 



4 A SIMPLE PARAMETRIZATION OF THE EXPLOSION 
PROPERTIES 

We next incorporate the properties of the explosion into our model. 
This can be done by relating the two quantities describing the ejecta 
profile, g and m, to the explosion parameters. In practice, obtaining 
such a relation requires a specific model for the fast material lead- 
ing the ejecta, and we use the standard model of |Matzner & McKee| 
( |1999) l. In this model, the original density in pre-explosion outer 
layers of the star, po, is a simple function of the relative distance 
from the edge of the star. 



Po(ro) = ; 5 = 1 - ro/R* 



(16) 



where i?* is the outer radius of the progenitor and and rg is the 
radial position within the progenitor. The power n is a constant 
which depends on the assumed physics of the envelope - n — 3 
for a convective envelope and n = 3/2 for a radiative one, and the 
density scale, p*, depends on the details of the progenitor structure. 

Given this initial density profile, [Matzner & McKee] jI999[ > 
have shown that after the passage of the shock the velocity of the 
material in the very outer envelope follows an approximate depen- 
dence on the energy of the explosion, E, and the total mass of the 
ejecta, M: 



v{ro) 



1/2 



47r 



-I3n 



(17) 



where /3 « 0.19 is practically independent of the details of the pro- 
genitor and the explosion, and fp = p*/p, where p = SM/AtvRI 
is the average mass density within the star. The asymptotic value 
of the coefficient is ~ 2Aj,,s, where A^^s ~ 0.79 is the ap- 
propriate value for the shock velocity profile at the onset of the ex- 
plosion. If the wind were dynamically unimportant and the ejecta 
could reach reach this asymptotic velocity profile, then the resulting 
density profile obtains the form ( [Rabinak & Waxman|2010^ , 



p{5m,r) = 



M (n + 1) ^ 



(18) 



where 5m{S) is the Lagrangian mass fraction of an element of the 
ejecta which was initially at 5 in the pre-explosion profile. With 
5m = 3 f p5^"'^^^ / (n + 1) and approximating that an element with 
5m reaches a radius r at a time of t ~ r/n we arrive at an expres- 
sion of the ejecta in the form of equation (j2]l: 



p{r,t) = 



47r 



l/r. 



{A-.fE^l^- 



(19) 
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where 9 ^{n+ l)/(2/3n), and 



n + 1 



(20) 



Again, we assume that the ejecta is cold and has developed the 
density profile of equation l|2j. Joining equations \\9) and (|9]l sets a 
quantitative constraint on the combination of the explosion energy, 
E, and the ejecta mass, M. Specifically for /? = 0.19 and n = 3 
we have m ~ 10 and 6 ~ 3.5, and with the aid of the characteris- 
tic values derived in the self-similar solution we aiTive (after some 
algebra) at the final result: 



E 



10^1 ergs 



M 



-5/7 



n 8 4-2 .2/21. -4f .-2/7 10/21 pl6/21 -20/21 

Although this equation is approximate, it implies that within the 
framework of our model, the plausible range of values for the en- 
tire prefactor in equation (j2TJ is rather limited. The values of x and 
Xe are tightly constrained in the self-similar solution, and do not 
introduce a significant uncertainty. The coefficient A f which does 
depend on the power of the ejecta density profile is suppressed by 
the 2/7 power, and the unknown factor fp is suppressed by the 2/21 
power Finally, ~ 1.6 for ejecta that has approached its asymp- 
totic velocity. The resulting equation can be applied as a sanity 
check for the value of the product EAI^^^'' , when a flash is sus- 
pected to be the result of a supernova-wind breakout scenario. 



5 THE X-RAY FLASH IN XRO080109/SN2008D 

We now examine the predictions of our model for the best can- 
didate to date of an X-ray flash due to shock breakout in a wind 
XRO080109, associated with the type Ibc supernova SN2008D 
( [Soderberg et al. 2008 ). Whereas a similar analysis can be applied 
to the case of GRB060218/SN2006aj, this source requires relativis- 
tic expansion ( |Li|2007[), and possibly a s ignificant deviation from 
spherical symmetry i Waxman et al.|2007] ), so we do not consider it 
here. 

The bright X-ray transient XRO080109 was serendipitously 
discovered during a scheduled Swift Telescope observation of the 
galaxy NGC2770 (at a distance of ~ 27 Mpc). The observed spec- 
trum shows a non-thermal shape with a power-law frequency de- 
pendence of the photon number flux per unit frequency, N^, cc i/"^ 
with r = 2.3 ± 0.3, through the observed range of X-ray energies, 
0.3-lOkeV. The follow-up observation by the Ultraviolet/Optical 
Telescope (UVOT) on board Swift discovered a counterpart Type 
lb supernova, denoted SN2008D. Analyses of the main lightcurve 
of SN2008D ( |Mazzali et a l. 2008 ; M odjaz et al. 2009 , Tan aka et al.| 
|2009[ > favor an underlying explosion of a compact progenitor, pre- 
sumably a Wolf-Rayet star due to the lack of hydrogen lines in the 
lightcurve. A consistent result was found by Rabin ak & Waxman| 
( |2010) in their analysis of the early UV/optical lightcurve, in which 
they inferred a initial progenitor radius of ~ lO^^cm. Note that for 
such a compact progenitor, a typical wind velocity should be of 
order 1000 km s"^ (similar to the escape velocity). 

The lightcurve of XRO080109 displayed a rapid rise and an 
exponential decline, with an e-folding time in the declining phase 
of approximately te = 130s. Such a timescale is too long to be 
consistent with breakout through the surface of a bare star, and the 
most plausible explanation links this decline with the diffusion of 
the radiation from the shock through the unshocked wind, te — 



tdis, which we will use below. For an alternative interpretation, that 
the lightcurve was shaped by an aspherical explosion, see |Couch| 
|etaL| i |20T0) . 

Quantitative estimates of the total energy of the burst (assum- 
ing isotropic emission) depend on the modeling of the column den- 
sity near the source and extinction; |Soderberg et al.| ( |2008^ esti- 
mated the total energy in the burst to be Ex ~ 2 x 10*''ergs, for 
which the appropriate breakout radius arising from equation iTTSj is 
7?i30 ~ 1.6 X lO^^cm. For the lower value of £;x ~ 6 x 10*^rgs 
inferred by [Modjaz et al.| ( |2009) , the breakout radius is reduced 
to Rbo ~ 1-1 X 10 cm. In either case, the breakout must have 
occurred in a wind rather than on the surface of a bare Wolf- 
Rayet progenitor. Our result is consistent with t he lower limit of 



7 X 10 cm for the breakout radius found by Soderberg et al. 



( |2008^ , based on the fact that the thermal component of the X-ray 
flash lies below the XRT low energy cutoff of ~ 0.1 keV. 

We note that the diffusion time scale indicates a density pro- 
file with K ~ 1.95 X lO^^g cm"^, corresponding to AI/v^ ~ 
3.85 X lO-^MQyr^VClOOOkms"^). This value is an order of 
magnitude larger than the range estimated by [Soderb erg et al.| 
( 2008 1 and Che valier & Fransson| j2008| > from the radio observa- 
tions of SN2008D, assuming a 10% efficiency in magnetic field 
production behind the shock front. Taken at face value, our result 
implies a lower efficiency for magnetic field production, or an en- 
han ced mass loss rate just p rior to explosion, or both. We note that 
that Soderberg et al. 



2008 1 estimated a low A4 also from the break- 



out flash, but this discrepancy is alleviated once our value for the 
breakout radius is used (rather than their lower limit), along with 
with a specific opacity of k — 0.2 cm^g^^ and a remaining optical 
depth of r =2-3 in the unshocked wind. 

The corresponding total amount of mass involved in the for- 
ward and reverse shocks at breakout comes out to be of order 
IQ-^A/q. Although small, this mass is several orders of magnitude 
larger than that enclosed in the last fe w optical depths of a bar e 
Wolf-Rayet star with ^ lO^^cm ( [Matzner & McKee||l999"l ), 
again supporting the wind breakout scenario. 

Finally, the combination of explosion energy and ejected mass 
in equation ^ yields (£/10^^ergs)(M/MQ)-^/'^ ^ 0.26 - 
-0.64/p/^^ for the total energies of Ex = 0.6-2 x lO'^^ergs. 
Depending on the value of fp, this result appears to lie between 
the combination originally suggested by .Soderberg et al., ^2008| 
of (£/10^^ergs) = 2, (M/Mq) = 5, and the more energetic 
explos ions inferred later by |Mazzali et al.| ( |2008| ) and [Tanaka et al.| 
( 2009| , (£;/10^^ergs) = 6, (M/Mq) = 6-8. 

We caution that {Rbo/R*) ^ ll-16(i?v,/10"cm)-^ is 
modest, and so the deviation from the cold ejecta model may be 
significant (and even more so if the progenitor radius was a few 
lO^^cm). This deviation does not change much the energetics of 
the model, which is why the general results are consistent with 
shock breakout in a wind, but the details of the shock structure are 
most likely different than those predicted by our solution. A further 
complication will arise from the wind density profile, which close 
to the star may be significantly steeper than up ~ r^"^ depen- 
dence ( [Li|2007^ . If the shocked wind material is dominated by this 
component, additional quantitative corrections must be included in 
our model. Such a wind structure may also provide a natural ex- 
planation the relatively large mass loss rate we required assuming 
p ^ throughout the entire wind. 

Another point of interest is that the shock velocity at break- 
out we find for SN2008D is quite large, with /3 ^ 0.3-0.4. This 
range is still consistent with our non-relativistic approach, and is 
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also in agreement with the velocity range inferred from later VLBI 
measurements ^Bietenholz et al.||2009[ l. However, given the rela- 
tively small breakout radius, it is most likely an upper limit: at 
{Rbo / R*) ~ 10-15 the ejecta may still be carrying a signifi- 
cant fraction of its initial internal energy that was not converted to 
kinetic energy. We also emphasize that if the optical depth at break- 
out is only r — j3^^ ~ 2.5 — 3, the diffusion approximation be- 
comes questionable. Finally, at /3 > 0.2 pair creation in the shock is 
non-neghgible ( Weaver] 1976[|Katz et al.|2010[|Budnik et al.|2010| l, 
causing the equation of state to deviate from a pure 7 = 4/3 ideal 
gas. Nonetheless, our simple model is robust enough to demon- 
strate that the shock velocity in SN2008D was large enough for the 
photon energy distribution not to be Planckian, leading naturally to 
the observed non-thermal X-ray spectrum. 



6 CONCLUSIONS AND DISCUSSION 

We have presented a simple model for the properties of shock 
breakout in a supernova explosion embedded within an optically 
thick wind. By assuming that the ejecta and wind form a forward- 
reverse shock structure ( |Chevalier|1982^ , we have shown that the 
breakout radius can be estimated based on the observationally de- 
terminable values of the total energy in the flash and the diffusion 
time scale from the shock to the wind's photosphere. Among other 
things, we demonstrated that the time scale for photon diffusion in 
the shocked region is shorter than the time scale for photon diffu- 
sion from the shock to the photosphere, justifying the underlying 
assumption that breakout will manifest itself as a brief flash even 
in the presence of an optically thick wind. We also found that if 
the ejecta time-dependent density profile is derived from the initial 
progenitor model of |Matzner & McKeeH1999) , the combination of 
explosion energy and ejecta mass can be estimated from the ob- 
served quantities of the shock breakout. For a low flash energy and 
long diffusion time, a blackbody spectrum is expected to emerge. 
In the opposite regime, free-free emission cannot generate a suf- 
ficient number photons to establish a thermal spectrum during the 
shock passage (Katz et al. 2010 ). 

While limited by its simplifying assumptions, our model de- 
scribes the general features of a shock breakout in a wind around 
an exploding Wolf-Rayet star, characterizing a Type Ibc supernova. 
We applied our model to the X-ray flash observed in association 
with SN2008D ( So defberg et aH2008) l and found sensible values 
for the shock breakout radius (~ 1.1-1.6 x lO^^cm) as well as 
the properties of the entire explosion. However, such a breakout ra- 
dius is probably not sufficiently larger than the progenitor radius 
for our model to be fully applicable. It is therefore possible that we 
are overestimating the shock velocity at breakout, which we found 
to be ~ 0.3-0.4. Nonetheless, our model clearly demonstrates 
that if XRO080109 was indeed shock breakout in SN2008D, the 
shock velocity was certainly large enough to prevent a blackbody 
spectrum from developing, in agreement with observations. 

Additional note: After the completion of this paper a related 
preprint was posted by |Chevalier & Irwin| ( |201 l"| l. They too use the 
self-similar solution to assess the total energy available for shock 
breakout through a wind, but their work focuses on a different set- 
ting involving a red supergiant exploding into a much denser wind. 
In their case the typical breakout radius is Rbo ~ lO^^cm and the 
mass of the ejecta involved in the shock comes out to be > 1Mq\ 
hence the assumption of a simple power law dependence for the 
density of the ejecta (equation|2]l may not be adequate (see §|3]l. 
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